The effect of uniaxial pressure on the magnetic anisotropy of the Mni2-Ac single-molecule magnet 
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We study the effect of uniaxial pressure on the magnetic hysteresis loops of the single-molecule magnet Mni2- 
Ac. We find that the application of pressure along the easy axis increases the fields at which quantum tunneling 
of magnetization occurs. The observations are attributed to an increase in the molecule's magnetic anisotropy 
constant D of 0.142(l)%/kbar. The increase in D produces a small, but measurable increase in the effective 
energy barrier for magnetization reversal. Density-functional theory calculations also predict an increase in the 
barrier with applied pressure. 
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I. INTRODUCTION 

Single-molecule magnets (SMMs) are fascinating systems 
in which each molecule behaves as a single rigid, high- 
spin object. Most SMMs have a large magnetomolecular 
anisotropy that impels the spin to point along a particular axis, 
the so-called easy axis. Thus, the spin's energy landscape can 
be described by a double- well potential, with an energy barrier 
between "up" and "down" orientations, as shown in Fig.[T] 

Much of SMMs' behavior can be understood from the sim- 
ple Hamiltonian^ 
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Here the first and second terms create the double- well poten- 
tial. The spin vector minimizes its energy by pointing along 
or antiparallel to the easy z axis. The constants D and A give 
the strength of the anisotropy, producing a classical barrier of 
height U = DS^ + AS"^ at zero field. The third term is pro- 
duced by the Zeeman interaction, with g = 1.96. When a 
field is applied along the z axis, the potential is tilted, making, 
say, the "up" direction lower in energy than the "down" direc- 
tion. W contains terms that do not commute with Sz and can 
therefore induce tunneling. 

Because SMMs have a finite spin S (~10), there are 25 + 1 
magnetic orientation states, m = —5, —S + 1, . . . , 5, associ- 
ated with the spin (levels in Fig.[T]). As the field along the easy 
axis is increased, levels in one well move down while levels 
in the other well move up. At certain fields^. 
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where N = —(m ^ m') = 0, ±1, ±2, . . ., levels m and m' in 
opposite wells align, allowing the spin to tunnel between wells 
and consequently producing a dramatic increase in the rate at 
which the spins come to equilibrium^. This effect, which is 
manifested by the appearance of steps in the hysteresis loops 
of most SMMs, has now been seen in hundreds of molecules 
and some rare-earth ions^ 

Mni2 acetate (henceforth referred to as Mni2) is the 
most studied SMM. Crystals of Mni2 typically contain two 



species^: a fast relaxing minority species and a slower re- 
laxing majority species. Here we exclusively study the lat- 
ter, which has a spin S — 10, a large energy barrier of ~70 
K, and is well described by Eq. [T] with D = 0.556 K and 
A = 1.13 mK^. Each magnetic molecule in a crystal is iso- 
lated from its neighbors, resulting in a total magnetic moment 
that is, to good approximation, an ensemble average of that of 
a single Mni2 molecule. At temperatures above ~ 0.5 K, the 
relaxation between the two wells is thermally assisted, with 
tunneling taking place from excited levels that are thermally 
populated. In this regime, the magnetic relaxation rate T is 
well described by an Arrhenius law^: 
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where Uef / is the effective barrier, defined as the energy dif- 
ference between the lowest level in the metastable well and the 
level m where the preponderance of tunneling takes place, i.e. 
Ueff = D{S^- m^j+A {S^ - m^)^gfiol^BHz{-S-my, 
ujQ, the Arrhenius pref actor, is a constant of order 10^ rad/s. 




m = 10 



FIG. 1 . Double- well potential diagram illustrating the energy levels 
for an S = 10 system. A magnetic field tilts the potential. When the 
field brings levels in opposite wells into alignment, the spins reverse 
via a process of resonant tunneling. At temperatures above ^-^0.5 K, 
this process is thermally assisted and tunneling occurs between levels 
near the top of the barrier, as indicated by the arrow. 

There have been numerous hydrostatic pressure experi- 
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ments on Mni2 and related molecules . Inelastic neutron 
scattering experiments have shown a few-percent increase in 
D when pressure of the order of a dozen kbar is applied^. In 
contrast, Levchenko et ai^ reported that pressure increases the 
measured magnetic relaxation rate, consistent with a decrease 
in D. In addition, studies have found that pressure converts 
some of the molecules between the fast and slower relaxing 
varieties, although the precise behavior of this change remains 
unclear^'^'^^'^^. In the present work, we apply uniaxial pres- 
sure along the easy axis of a single crystal of Mni2 and ex- 
amine the energy-level structure through measurements of the 
sample's hysteresis loops. We find that the pressure alters the 
field at which tunneling occurs, an effect we interpret as due 
to an increase in D. We also find a slight decrease in F with 
pressure, which we attribute to the increase in Uef / brought 
about by the increase in D. Our technique allows us to vary 
the pressure in situ at low temperature so that the sample is 
never under high pressure at room temperature, minimizing 
any structural distortions and deformations that may occur at 
such temperatures. 



II. EXPERIMENTAL SETUP 

Our uniaxial pressure apparatus is based on a design by 
Campos et alJ^. A diagram of the low-temperature portion of 
the apparatus is shown in Fig.O Pressure is applied through a 
pneumatic piston connected to a nitrogen gas cylinder, and is 
transferred to the sample cell at the end of the probe by non- 
magnetic stainless steel rods, capped by a small aluminum 
plug. This plug and the bottom cap of the cell press on two 
GIO phenolic "fingers" that transmit the pressure to a single 
crystal of Mni2, which is encased in a small (non-circular) 
cylindrical piece of epoxy (Stycast 1266), as shown in Fig.O 

The sample was prepared by placing it in freshly mixed 
epoxy within a Teflon mold. A 3-Tesla field was applied to 
align the crystal's easy axis (crystallographic c axis) with the 
axis of the cryostat's field. Once cured, the epoxy was ma- 
chined to make the top and bottom horizontal faces flat and to 
create a flat vertical face with the sample close to the surface, 
allowing coupling to a Hall bar detector. 

The bottom portion of Fig. [2] illustrates the bracket that 
holds the Hall sensor chip and its mounting board, the GIO 
fingers, and the epoxy "pellet" containing the sample. When 
everything was in place, one edge of the sample was aligned 
with one Hall sensor; another sensor relatively far from the 
sample was used to measure background signals. The signals 
from the two Hall bars were subtracted using an analog cir- 
cuit in order to cancel out common background signals. A 
regulated thermometer was placed close to the epoxy-sample 
pellet for temperature measurement and control. 



III. RESULTS AND ANALYSIS 

We acquired magnetic hysteresis loops from the sample at 
temperatures between 1.8 K and 2.2 K in 50 mK increments 
and with applied pressures, P, as high as ~1.8 kbar. For all 




FIG. 2. (top) Cross section schematic of our uniaxial-pressure cell, 
(bottom) Exploded view of the cell core, showing Hall detector chip 
mounted on circuit board, the bracket, GIG fingers, and the epoxy 
pellet containing the embedded sample. 



loops, the field sweep rate was dH/dt = 3.0 mT/s. A typical 
hysteresis loop (taken at 2.0 K and with no applied pressure) 
is shown in the inset of Fig.O The magnetization steps char- 
acteristic of resonant tunneling are apparent. The loop also 
displays features near zero field, which are the paramagnetic 
response of the fast-relaxing species as well as of some mag- 
netic impurities in the apparatus' materials. The main figure 
shows a close-up of the A/^ = 4 step for this loop. For every 
tunneling step clearly visible in each loop, we empirically fit 
the data with a hyperbolic tangent plus a straight line (to ad- 
just for residual background effects) to determine Hn{T^ P), 
the field center for step N (the field at which the tunneling rate 
is maximum). The main part of Fig. [3] shows an example of 
such a fit, superimposed on the data. For a given N and T, the 
step centers generally shift toward higher field with increas- 
ing pressure. Figure |4] shows Hn as a function of pressure for 
the = 4 step at 2.0 K. With the exception of cases where 
the step is very small, with a correspondingly small signal-to- 
noise ratio, Hn appears to depend linearly on pressure over 
the range of pressures studied. The results show no apparent 
baric hysteresis as evidenced by the fact that some of the data 
at ~0.45 kbar and all of the data at and ~0.67 kbar were 
taken after applying the highest pressure. 

The values of Hn depend on both pressure and tempera- 
ture. The latter dependence, primarily due to the effects of 
mean dipolar fields is significantly stronger than the effect 
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FIG. 3. Detail of the N = 4 step with fit. The data was taken at 2.0 
K and with no pressure appUed. The inset shows the full hysteresis 
loop from which the detail was taken. 



of pressure. To determine the effect of pressure on for a 
particular step, we simultaneously fit all of the data for that 
step. We take Hn for each step to be a linear function of 
pressure with a different zero-pressure intercept for each tem- 
perature: 



Hn{T,P) = Hn{T,0) 



N 



dP 



(4) 



Fitting all data for each step to Eq. (H we obtain a value 
of OHn /dP for that step. Figure |5] shows the values of 
OHn/OP as a function of N. (Here and below we have 
grouped together data for N and —TV.) The figure also shows 
a fit to a line with zero intercept. This result suggests that 
the anisotropy parameter D is increasing linearly with applied 
pressure, consistent with Eq. |2l which predicts dH^ / dP = 
-dD/dP when dA/dP = 0. From the fit, we obtain 
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a pressure-induced relative change of D of +0. 142(1 )%/kbar. 
This value differs somewhat from those found in experiments 
employing hydrostatic pressure^'^^. The difference may be 
due to the uniaxial nature of the applied pressure in our study. 

We determined the values of m (and m' = N — m) for each 
step by carefully ascertaining the values of Hn at zero pres- 
sure and, after correcting for dipole-field effects employing 
Eq.[2l From this analysis we find that the pair of levels (m, m') 
between which tunneling occurs are three levels below the top 
of the classical barrier for all three steps {N = 3,4,5) that we 
measured^^. 

As a consistency check of the above interpretation that pres- 
sure is inducing a change in I), we looked at how the pressure 
changes the effective barrier height Ueff- To do this, we also 
measured the maximum slope dM/dH for each step and the 
corresponding value of M at that point in the step. Under the 
approximation that the system has a single relaxation rate F, 



CD 

O 



18.08 



18.07 



18.06 



18.05 




18.04 



0.5 1 

P (kbars) 

FIG. 4. Step center field values as a function of pressure for the N : 
4 step at 2.0 K. The line is the result of a linear fit of the data. 
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FIG. 5. OHn I dP as a function of step number. The line is the result 
of a linear fit with fixed zero intercept. 



the relaxation is governed by the simple differential equation 

dM dM dH 



dt 



dH dt 



-T{M-Meq). 



(5) 



where the field sweep rate dH/dt is constant for our exper- 
iments. Meq is the equilibrium magnetization and is nearly 
equal to the saturation magnetization Mgat for the values of 
N and T investigated. With that approximation, F can be de- 
termined from Eq. [5] and the measured values of dM/ dH and 
M for each step: 
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At the same time, assuming a linear pressure dependence 
of Ueff (i.e. Ueff = Ueff,o + f/^/^P), Eq.|3]yields 



r = ro (iJ,T)e-^«//'P/^, 



(7) 



where the prime indicates differentiation with respect to P 
and To (i^, T) is the relaxation rate at P = 0. We neglect the 
small P dependence ujo may have. 

For each N and T, we determine F (P) using Eq. [6] and 
then fit the pressure dependence to Eq. |7]to determine t^g//- 
The results are shown in Fig. [6] as a function of T. The 
horizontal dashed lines are the expected values of U^^j^ = 
dD/dP [S'^ — m?) using the value of dD/dP determined 
from the results in Fig. [5] and the value of m we determined 
from the values of Hn for each N^^. 

The = 3 data is very scattered and we can only con- 
clude from it that the value of /7gjj agrees on average with 
the sign and order or magnitude of the prediction. For N = 4, 
the scatter is smaller (note the change in vertical scale) and 
the agreement with the prediction is better, although there is 
an unexpected decrease with increasing T. For N = 5, the 
data is consistently below the predicted value, but, again, of 
the same order of magnitude. The disagreements most likely 
arise from the assumptions of a single relaxation rate F (Eq.O 
and monodispersed values of D and A. It is well known that 
Mni2Ac has a distribution of relaxation rates due to a distri- 
bution in anisotropy parameters as well as tilts of the easy axis 
away from the crystallographic symmetry axis^^~-. These 
effects result in a broadening of the tunneling steps that ef- 
fectively creates a maximum value of dM/dH. Thus, when 
the steps are large (higher-T data for = 4 and all data for 
N = 5), the value of F estimated using Eq.[6]is nearly satu- 
rated, suppressing its P dependence. 

While it is hard to draw strong conclusions from the re- 
sults in Fig. [6l it is clear that this data is generally consonant 
with the interpretation that pressure is inducing a change in D 
and consistent with the hypothesis that dA/dP = 0. If we 
consider that A is also changing with pressure, it is difficult 
to devise a situation with significantly better agreement with 
the data. For example, allowing dA/dP = 0.5%/kbar (re- 
ducing dD/dP to 0.116(l)%/kbar) results in a very modest 
improvement in the fit in Fig. [S] but roughly doubles the pre- 
dicted values of /7gjj in Fig. [6l If we take an extreme case 
and fix dD/dP = 0, then a fit of the data in Fig. [5] results in 
dA/dP = 2.73(1) %/kbar, which predicts U'^ff to be in the 
range 0.27-0.29 K/kbar, clearly inconsistent with the results 
in Fig. [6l We cannot rule out that pressure is inducing small 
changes in A, but the only cases in which changes in A would 
produce significant effects require that its relative change be 
larger than the relative change in 1), a somewhat unlikely sce- 
nario. 

To complement the experimental findings, we performed ab 
initio calculations of the pressure dependence of the energy 
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FIG. 6. U'^ff as a function of T for N = 3, 4, 5, as labeled. The 
horizontal dashed lines are the expected values of U^ff using the 
value of dD / dP determined from the results in Fig. |5] Note that the 
three panels have different vertical scales. 



barrier of Mni2. We calculated the electronic structure and 
magnetic properties of the molecule using density-functional 
theory (DFT) including spin-orbit coupling (SOC). We used 
the plane-wave-based DFT code VASP^^-^^ and, for the 
exchange-correlation potential, we used the Perdew-Burke- 
Ernzerhopf generalized-gradient approximation^^. Projector- 
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augmented-wave pseudopotentials were used. The kinetic- 
energy cutoff was set to 580 eV. For our calculations, we used 
a unit cell of dimensions 24 x 24 x 24 containing one Mni2 
molecule. We performed self-consistent DFT calculations in- 
cluding SOC until the total energy converged to within 0.5 
/ieV. To simulate the effect of uniaxial pressure, we assumed 
uniform strain along the z axis of the molecule. This was im- 
plemented by contracting all z coordinates of the molecule by 
a given factor (up to 1.5%) while leaving the x and y coordi- 
nates unchanged. For each value of strain, we computed the 
magnetic anisotropy barrier U in zero magnetic field. We con- 
verted strain to applied pressure using the method discussed 
jj^27,28 (tantamount to a Young's modulus of 137.1 kbar). The 
results of these calculation are shown in Fig. [7] 
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FIG. 7. Calculated energy barrier as a function of uniaxial pressure. 
Calculations were performed using density-functional theory includ- 
ing spin-orbit coupling, as described in the main text. The line shows 
the result of a linear fit to the calculated values. 

These calculations indicate a qualitatively similar increase 
in U with pressure as observed experimentally. Fitting the 
data in Fig. [7] to a straight line, yields a slope of dU/dP = 



0.30(4) K/kbar, about a factor of four larger than predicted: 
dU/dP = S'^dD/dP = 0.0790(5) K/kbar. The discrepancy 
is not surprising: It indicates that most of the induced strain is 
intermolecular, which is explicitly not considered in the the- 
oretical model. In other words, the Young's modulus for the 
crystal is significantly smaller than that for the molecule itself. 



IV. SUMMARY 

We found that by applying uniaxial pressure along the easy 
axis of Mni2 acetate we were able to increase the magnitude 
of the applied field values at which resonant tunneling occurs. 
We explain this as an increase in the anisotropy parameter D 
in the system's spin Hamiltonian. We determine that pressure 
causes a change in D of 0.142(l)%/kbar. This interpretation 
is consistent with the pressure-induced changes in the relax- 
ation rate F as well as ab initio calculations of the the effect 
of pressure on the barrier height U. 

While in this experiment, the pressure was applied along 
the molecule's symmetry axis, one could also apply pressure 
along a perpendicular direction, which could affect the tun- 
neling rate by breaking the symmetry of the molecule. In fact, 
a theoretical calculation predicts that pressure applied along 
the hard axis of a four-fold symmetric molecule could modu- 
late the tunneling rate via a geometric-phase effect2^. Future 
experiments may attempt to observe this effect. 
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